Gravitational wave signatures of the absence of an event horizon. 
I. Nonradial oscillations of a thin-shell gravastar 
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Gravitational waves from compact objects provide information about their structure, probing 
deep into strong-gravity regions. Here we illustrate how the presence or absence of an event horizon 
can produce quahtative differences in the gravitational waves emitted by ultra-compact objects. In 
order to set up a straw-man ultra-compact object with no event horizon, but which is otherwise 
almost identical to a black hole, we consider a nonrotating thin-shell model inspired by Mazur 
and Mottola's gravastar, which has a Schwarzschild exterior, a de Sitter interior and an infinitely 
thin shell with finite tension separating the two regions. As viewed from the external space-time, 
the shell can be located arbitrarily close to the Schwarzschild radius, so a gravastar might seem 
indistinguishable from a black hole when tests are only performed on its external metric. We study 
the linearized dynamics of the system, and in particular the junction conditions connecting internal 
and external gravitational perturbations. As a first application of the formalism we compute polar 
and axial oscillation modes of a thin-shell gravastar. We show that the quasinormal mode spectrum 
is completely different from that of a black hole, even in the limit when the surface redshift becomes 
infinite. Polar QNMs depend on the equation of state of matter on the shell and can be used 
to distinguish between different gravastar models. Our calculations suggest that low-compactness 
gravastars could be unstable when the sound speed on the shell Vs/c > 0.92. 
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I. INTRODUCTION 

Black holes (BHs), once considered an exotic mathe- 
matical solution to Einstein's field equations, have now 
been widely accepted as astronomical objects P, d, H, 3- 
Stellar-mass BHs arc believed to be the final stage of 
the evolution of sufficiently massive stars. Massive BHs 
seem to populate the center of many galaxies at low red- 
shift, and must have played an important role in the 
formation of structure in the universe. Most evidence 
supporting the astrophysical reality of BHs comes from 
the weak-gravity region, i.e. from observations probing 
the space-time several Schwarzschild radii away from the 
event horizon. Attempts to rule out possible alterna- 
tives to BHs usually rely on general relativity being the 
correct theory of gravity, and/or on constraints on the 
equation of state of matter at high densities. For mas- 
sive BHs the most precise measurements so far come from 
observations of stellar proper motion at the center of our 
own galaxy, indicating the presence of a "dark object" of 
mass M ~ (4.1 ± 0.6) x 10*^ M© Recent millimeter 
and infrared observations of Sagittarius A* , the compact 
source of radio, infrared and X-ray emission at the center 
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of the Milky Way, infer an intrinsic diameter of 37_2o mi- 
croarcseconds, even smaller than the expected apparent 
size of the event horizon of the presumed BH Q . Some of 
the exotic alternatives to a BH (such as "fermion balls" ) 
are incompatible with the observations 0] and any dis- 
tribution of individual objects within such a small region 
(with the possible exception of dark matter particles or 
asteroids, which however should be ejected by three-body 
interactions with stars) would be gravitationally unstable 
@, Q. In a recent attempt to probe the event horizon, 
Broderick and Narayan have analyzed the observations 
of Rcf. If the object at the center of our galaxy had 
a surface it would be hot enough to glow with a steady 
emission of infrared light, but no such glow has been de- 
tected [lo| . This and similar arguments are inevitably 
dependent on the gas distribution and on details of the 
accretion process, and they really set lower limits on the 
gravitational redshift corresponding to the hypothetical 
surface replacing the event horizon (see e.g. [11| for a re- 
view). Indeed, some hold the view that an observational 
proof of the existence of event horizons based on electro- 
magnetic observations is fundamentally impossible p^ . 
For these reasons, model-independent tests of the strong- 
field dynamics of BHs and studies of the gravitational 
radiation signatures of event horizons are necessary to 
confirm or disprove the BH paradigm p^ . [l^ . 

Gravitational wave detectors offer a new way of ob- 
serving BHs, complementing the wealth of information 
from present electromagnetic observations As first 
proposed by Ryan [l^, , an exquisite map of the ex- 
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ternal space-time of BHs (outside the innermost stable 
orbit, if there is any) can be constructed by observing 
the gravitational waveform emitted when a small com- 
pact object spirals into the putative supermassive BH 
at the center of a galaxy with the Laser Interferometer 
Space Antenna (LISA). As an extension of Ryan's work, 
Li and Lovelace considered the small object's tidal cou- 
pling with the central object and showed that by study- 
ing details of radiation reaction, information about the 
space-time region within the orbit can also be obtained 
[l7| . Ryan's proposal to map space-times using inspi- 
ral waveforms is promising, but the data analysis task is 
affected by a "confusion problem" : the possibility of mis- 
interpreting truly nonKerr waveforms by Kerr waveforms 
with different orbital parameters [isl . fioj . This ambigu- 
ity was shown to be resolvable if the orbit is known to 
be circular or if one only probes the mass, spin and 
quadrupole moment of the object using waveforms gener- 
ated in the weak-gravity region [20| . A different approach 
to test the BH nature of ultra-compact objects is based 
on measuring several of their free oscillation frequencies 
( "ringdown waves" ) and comparing them with the quasi- 
normal mode (QNM) spectrum of BHs [2l|. These tests 
of the "no-hair theorem" require a signal-to-noise ratio 
which could be achieved by advanced Earth-based grav- 
itational wave interferometers and they are one of the 
most promising science goals of LISA [22|, [2^ [^J] ■ 

In this paper we explore how to test possible alterna- 
tives to the BH paradigm. We retain the "conservative" 
assumption that general relativity is the correct theory 
of gravity and we focus on possible tests of the existence 
(or absence) of an event horizon. 

Theorists conceived several families of compact objects 
with no event horizons. For example, boson stars are 
horizonless compact objects based on plausible models 
of particle physics at high densities, and they are (still) 
compatible with astrophysical observations [H, [2^ . Be- 
ing indistinguishable from BHs in the Newtonian regime, 
boson stars are good "strawnien" for supermassive BHs. 
The space-time of nonrotating spherical boson stars can 
approximate arbitrarily well a Schwarzschild geometry 
even close to the event horizon, and being very compact 
it is not easily distinguishable from a BH by electromag- 
netic observations [26|, . Building on Ryan's proposal, 
Kesden et al. showed that the inspiral of a small compact 
object into a nonrotating boson star will emit a rather dif- 
ferent gravitational waveform at the end of the evolution, 
when the small object falls into the central potential well 
of the boson star instead of disappearing into the event 
horizon of a BH [28J. Several authors have computed 
the QNM spectrum of boson stars, showing that it is re- 
markably different from the QNM spectrum of BHs and 
lending support to the feasibility of no-hair tests using 
QNM measurements [H, HO, IMi . 

Another proposed alternative to massive BHs, which 
we shall focus on in this paper, are the so-called gravas- 
tars [3^ . The gravastar model assumes that the space- 
time undergoes a quantum vacuum phase transition in 



the vicinity of the BH horizon. The model effectively 
replaces the BH event horizon by a transition layer (or 
shell) and the BH interior by a segment of de Sitter space 
[33I . |33 |. Mazur and Mottola argued for the thermo- 
dynamic stability of the model. A dynamical stability 
analysis by Visscr and Wiltshire [1^ confirmed that a 
simplified version of the gravastar model by Mazur and 
Mottola is also stable under radial perturbations for some 
physically reasonable equations of state for the transition 
layer. Chirenti and RezzoUa [3^ first considered nonra- 
dial perturbations of gravastars, restricting attention to 
the relatively simple case of oscillations with axial parity. 
They computed the dominant axial oscillation modes and 
found no instabilities. In analogy with previous studies 
of the oscillation modes of boson stars [2^, HO, [3l| , they 
confirmed that the axial QNM spectrum of gravastars 
can be used to discern a gravastar from a BH. In the 
thin-shell limit, the axial QNM frequencies of Ref. [sOl 
and our own calculations recover Fiziev's calculation of 
the axial QNMs of ultracompact objects with a totally 
reflecting surface [s^l- 

In this paper we study the stability with respect to 
nonradial oscillations of a simplified "thin-shell" gravas- 
tar, retaining most of the essential features of the original 
model. We consider both axial perturbations (reproduc- 
ing and extending the results of Ref. [11] ) and polar per- 
turbations. In the polar case, the matching of interior 
and exterior perturbations at the gravastar shell requires 
a more careful analysis because (unlike the axial case) 
polar perturbations of spherical objects actually induce 
motions of matter, which in turn couples back to gravi- 
tational perturbations. For this reason, nonspherical po- 
lar perturbations provide a more stringent test on the 
gravastar's overall stability. Polar perturbations are also 
crucial in studying the dynamics of objects orbiting the 
hypothetical gravastar. In order to treat polar pertur- 
bations of a nonrotating thin-shell gravastar we set up 
a rather generic formalism combining standard pertur- 
bation theory (in the Regge- Wheeler gauge) with Israel's 
junction conditions [ssj . The formalism can be applied to 
gravitational perturbations of any spherically symmetric 
space-time characterized by regions with different cosmo- 
logical constants separated by infinitely thin shells with 
finite surface energy and tension. Because matter is con- 
centrated on these shells, the junction conditions deduced 
here will be sufficient in describing the linear dynamics 
of matter. Quite predictably, these conditions depend on 
the equation of state of the shell. As an application of 
the formalism we study the QNM spectrum of polar and 
axial perturbations of gravastars, exploring the nonra- 
dial stability of these objects. Polar QNMs (unlike axial 
QNMs) depend on the equation of state of matter on the 
shell: they can be used not only to distinguish a gravas- 
tar from a BH, but also to distinguish between different 
gravastar models. We also find that the imaginary part 
of some QNMs seems to have a zero crossing when the 
gravastar is not very compact and the speed of sound on 
the shell is superluminal, suggesting that some gravastar 
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models may be unstable under nonradial perturbations. 
In a companion paper we will apply the formalism devel- 
oped in this paper to study gravitational radiation from 
compact objects inspiraling into nonrotating gravastars. 

The plan of the paper is as follows. In Section |TT] 
we review our static thin-shell gravastar model. Section 
mil sketches the calculation of axial and polar gravita- 
tional perturbations and of the matching conditions at 
the gravastar shell. Details of the matching procedure 
are provided in Appendix [X] and details of the QNM 
calculation are given in Appendix [Bl Our numerical re- 
sults for the polar and axial QNM spectra are presented 
in Section HVl and supported by analytical calculations in 
the high-compactness limit in Appendix [Cl 

We use geometrical units {G = c = 1). The Fourier 
transform of the perturbation variables is performed by 
assuming a time dependence of the form e"^*. Greek 
indices {fi , u , . . .) refer to the four-dimensional space- 
time metric. Latin indices (i ,j , . . .) refer to the three- 
dimensional space-time metric on the shell. Latin indices 
at the beginning of the alphabet (a , 6 , . . . ) refer to the 
spatial metric on a two-sphere. 



II. EQUILIBRIUM MODEL 

The metric for a static thin-shell gravastar has the form 

dsl = -f{r)dt^ + -T^dr^ + r\de^ + sin^ 6*^) . (2.1) 
h{r) 



Here 



r 



r > a. 



a h{r) = a ( 1 3^''^ ) ' ^ < a , 



(2.2) 



where M is the gravastar mass measured by an outside 
observer, and p = 3M/(47ra^) is the "energy density" of 
the interior region. The coordinate system (t, r, 6', ip) has 
been chosen in such a way that the thin shell occupies a 
coordinate sphere with r ~ a. The space-time reduces to 
de Sitter for r < a, and to Schwarzschild for r > a. The 
junction conditions on the r = a surface have already 
partially been chosen by requiring the induced metric to 
be continuous across the shell, which also dictates that 
/(r) be continuous at r = a, or 



, 2M ( ^Trpa^ 

1 = Q' 1 — 

a \ 6 



(2.3) 



In this paper we shall usually drop the dependence of /(r) 
and h{r) on r. From the jump in the radial derivatives 
of / we could easily obtain the two defining properties of 
this gravastar model: the surface energy density E and 
surface tension 0. The junction conditions read [3^: 



Si' 



S' 

7.. 2 



(2.4) 



where the symbol "[[ ...]]" gives the "jump" in a given 
quantity across the spherical shell (or r = a), i.e. 



[[A]]^A{a+)-A{a^). 



(2.5) 



The indices i and j correspond to coordinates i, 9, and 
Lp which parameterize curves tangential to the spheri- 
cal shell, Kij 

(0,1,0,0)/V5 

surface stress-energy tensor 



—ViUj is the extrinsic curvature, = 
is the unit normal vector, and Sij is the 



Sij = (S - Q)uiUj - Qjij , 



(2.6) 



where = ^^-1/^(4(1,0,0,0) (or u = yZ-l/gt^dt) is 
the four- velocity of mass elements on the shell and ^af3 = 
9ap — n-aiT-p is the induced 3-metric on the shell. We then 
have 



5 

= (s-e)i 



In the static, spherically symmetric case. 



~9ij,r 



(2.7) 



(2.8) 



Discontinuities in the metric coefficients are then related 
to the surface energy and surface tension as [ssj 



[[Vh]] = -47raE , 



f 



87r(S]-2e). (2.9) 



In order to summarize the above relations and reveal the 
independent parameter space of a thin-shell gravastar we 
define 



Ms = Aira^T, , 



(2.10) 



which would be the volume- and surface-energy contents 
of the gravastar. In terms of M^, Mg and a, we can 
obviously solve for p, E, and in addition we have 



M 



e = 



2M„ 



Ml 
2a 



My+Ms\ll 

1 - 2M/a 
1 - 2My/a ' 

1 r l-4M„/a 
Stto ^1 - 2My/a y/1 - 2M/a 



1 - M/a 



(2.11) 
(2.12) 

.(2.13) 



As a consequence, gravastar types can be specified by 
the dimensionless parameters My /a and Ms /a. In this 
paper we only consider a simplified version of the orig- 
inal Mazur-Mottola gravastar, which has vanishing sur- 
face energy (S = 0) and 



M = My, SnaO = - 



3M/a 



y/1 - 2M/a 



(2.14) 



However it is convenient to keep our notation general 
enough, because nonradial oscillations of a gravastar will 
in general produce nonzero variations of the surface en- 
ergy (i.e., 7^ 0). 
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III. GRAVITATIONAL PERTURBATIONS with 

111 both the interior (de Sitter) and exterior 
(Schwarzschild) background space-times we consider per- 
turbations in the Rcgge- Wheeler gauge [1^, writing 

ds^ = dsl + {Snwgtii^)dx^dx'^ (3.1) 



\\SRW9^L,^\\ = 



~'^a[t,r)^—— — ho[t,r)sm 



h{r) 



sin 6 dip 
1 dYi„, 



sin dip 



de 

dYun 

de 



/ii(i,r)sin( 


r2 sin^ eK{t,r)Yir, 



(3.2) 



where Yim{e,4>) denotes the ordinary spherical harmon- 
ics and "*" stands for terms obtainable by symmetry. In 
this gauge the perturbations split into two independent 
sets: the metric functions ho and hi are axial or odd par- 
ity perturbations, while Hq , Hi , H2 , K are polar or even 
parity perturbations. The linearized Einstein equations 
automatically require TJq = H2 ^ H . 

In the rest of this section we work out perturbations of 
the gravastar space-time, including the dynamics of the 
shell itself, in three steps. In Sec. IIII Al we present the 
well-known solution of the perturbation equations in de 
Sitter space-time in terms of hypergcometric functions 
and choose the solution that is regular at the origin [r = 
0). In Sec. IIII B I we review metric perturbations in the 
Schwarzschild exterior. Finally, in Sec. IIII Cl we work out 
the junction conditions relating the interior and exterior 
Regge- Wheeler perturbations. 



In the de Sitter interior both axial and polar perturbations 
can be reduced to the study of the (frequency-domain) 
master equation 



2,1/ in 



drl 



+ [uj'--VUr)\ 



0, 



r < a, 



where 



Mn(.)^^./(.) 



(3.4) 



(3.5) 



and wc introduced the tortoise coordinate, defined as 
usual by the condition dr/dr^, = /(?"), which in this case 
yields 



87rp 



-arctanh 



V 3 



r < a . 



(3.6) 



A. The de Sitter interior 

The usual way to obtain the interior solution for per- 
turbed stars is by direct integration of the system of 
ODEs [13, Sll, [121. Integrating a regular solution from 
the center would give boundary conditions at the stel- 
lar radius (or, in the case of a gravastar, at the location 
of the shell). For the de Sitter interior (r < a) no nu- 
merical integrations are required, because a regular so- 
lution of the perturbation equations can be obtained in 
terms of hypergcometric functions [1, \^ . To establish 
notation we review the basic equations below. Let us 
express the metric in terms of a compactness parame- 
ter C = {2M/a)^, related to the parameter /x = M/a of 
Ref. [Hi by C = S^i^. Then we have (assuming E = 0) 



In terms of r, the master equation reads 



Qj.2 



I dr 



P 



= , 



(3.7) 



where a prime denotes a derivative with respect to r. 
Near the origin, solutions of Eq. p.7p behave as Ar'+^ + 
Br^K By requiring regularity at the center (r = 0) we 
get, up to a multiplicative constant, 



r'+'{l-C{r/2MY)'"v^ 
l + 2-i^ l + l --2^^' 



2 ' 4M2 



,(3.8) 



where F{a,b,c, z) is the hypergcometric function [4J]. 
From 5*'" and its derivative we can get the axial per- 
turbation variables in the frequency domain [4^ : 



fir) = 1 - — r = 1 



2M 



l-C{r/2Mf. (3.3) 



hi 



i d 
LO dr^ 



(3.9) 
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The polar metric functions K and Hi can be obtained 
immediately from: 



R = -^i; + — 

2r dr^ 

iiur d^"' 

f \ r dr^, 



(3.10) 
(3.11) 



The quantity H{= Hq = H2) and its derivatives can then 
be found from the algebraic relation 



2 / / 

{l-l){l + 2) 



.l{l + l)Cr 



(3.12) 



This procedure fixes all metric quantities and their 
derivatives in the interior. Here and henceforth in the 
paper we drop the dependence of Hq, hi, H, Hi and K 
on uj and r. 



B. The Schwarzschild exterior 

In the Schwarzschild exterior, axial a nd p olar pertur- 
bations obey different master equations [4^. The deter- 
mination of the axial perturbation variables can still be 
reduced to the solution of the Regge- Wheeler equation 
[111, a Schrodinger-like ODE identical to Eq. ([XT)) : 



2iTiOUt 



52* 



Qj.2 

where 



/ dr 



^out('') ^out 



0, (3.13) 



(3.14) 
, with 



The metric can then be obtained from Eqs 
/(r) given by Eq. 

The perturbed Einstein equations relate the polar vari 
ables {K, H, Hi) via three differential equations: 



— ifHi)-iLu{H + K) = 0,(3.15) 
dr 



-iuoHi + f{H' - K') + J'H = 0,(3.16) 



K' - 



H 



L 
2/ 



^\-^Hi . 0,(3.17) 



and an algebraic relation 



1 



2 



1 



r_r_ 



.1(1 + 1] 

4uj 



f 



l{l 



2 2 



Hi 

H = 0. 



K 



(3.18) 



Note that if we make the appropriate choice for /(r). 
Eqs. (|3.15p - (|3.17p and the algebraic relation also apply 
to the interior de-Sitter spacetime [cf. Eq. (|3.12p ]. 

The Zcrilli function (r) [i^ , which satisfies a wave 
equation, and its spatial derivative are also constructed 
from Hi and K as 



^out 

dr^ 



Hf' 



A2 - A1A3 
A2 - A1A3 



where 



Ai 
A3 



GAP + A/2 (1 + A/2) 7-2 + 3/2Ai\/r 

r2 (3M + rA/2) 
iuj {-3AP - 3/2XAIr + 



r{UI + rX/2)f 



r 



(3.19) 
(3.20) 

(3.21) 

(3.22) 
(3.23) 



and A = (^ — l)(/ + 2). The metric perturbations can then 
be obtained by integrating the Zerilli equation outwards, 
starting from r = a+. It is also useful to recall that 
(in the exterior) we can switch from the Zerilli function 
Z™*(r) (|3:20l) to the Regge- Wheeler function ^'°"'(r) by 
using a differential relation between polar and axial vari- 
ables discussed in Chandrasekhar's book [491: 



where 



a± 



r]Z° 



A(A + 2) 
12M 



dZ° 



dr^ 



dr^ 



1] 



A(A 



12M 



6M/(r) 
r(Ar + 6AI) 



From and dZ°^'^ /dr^, we can easily compute and 
d'I'°"*/o?r^ outside the shell and use them as initial condi- 
tions to integrate the Regge- Wheeler equation outwards. 
Leins et al. [431 showed that this procedure is convenient 
to compute polar oscillation modes by the continued frac- 
tion method; more details on the QNM calculation are 
given in Appendix [B] 



C. The matching conditions 

In this section we discuss the most delicate part of 
the perturbation problem, namely the junction condi- 
tions for the Regge- Wheeler perturbation variables across 
the shell. Here we only outline our strategy and present 
the results; more details are given in Appendix \K\ Tech- 
nically, the application of Israel's junction conditions is 
easier if the shell's world tube happens to coincide with a 
fixed coordinate sphere at constant radius. However this 
is incompatible with choosing the Regge- Wheeler gauge 
in both the interior and the exterior, which is convenient 
to cast the perturbation equations into simple forms. In 
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fact, such a choice of gauge does not leave any freedom. 
We must exphcitly parametrize the 3-dimensional mo- 
tion of each mass element on the shell and then per- 
form the matching on a moving shell. In order to take 
advantage of both the simplicity of the field equations 
and the convenience of matching on a fixed coordinate 
sphere, we carry out the matching in the following way. 
We first construct a particular coordinate transformation 
(for both the exterior and interior space-times) such that 
in the new coordinate system, any mass on the shell re- 
mains static on the coordinate sphere with r = a. In 
this new coordinate system the metric perturbations will 
no longer be Regge- Wheeler, but will be augmented by 
quantities that carry information about how masses on 
the shell move in the Regge- Wheeler gauge. The stress- 
energy tensor of masses on the shell will correspondingly 
be modified. We then carry out the matching at r = a 
and obtain junction conditions relating the interior and 
exterior metric perturbations, plus equations of motion 
for matter on the shell. As could be anticipated from 
the general features of oscillations of nonrotating stars, 
axial perturbations do not couple to matter motion and 
the axial junction conditions are very simple, basically 
imposing continuity of the master variable 5* and of its 
first derivative. Polar perturbations, on the other hand, 
do couple to matter motion, so polar junction conditions 
do involve the shell dynamics, i.e. its equation of state. 

We parametrize the world line of matter elements on 
the shell in terms of the coordinates (t, r, 6, (p) as follows: 



t = 
r = 

e = 



(3.24) 



where and (p^ identify physical mass elements on the 
sphere, while r parametrizes their proper time. Note 
that the Lagrangian equations of motion will not be the 
same for the interior and exterior space-times. Therefore, 
points with the same t , 9 and ip coordinates are not in 
general the same when viewed from the interior and from 
the exterior. As shown in Appendix]^ the four- velocity 
of the mass element (0*,(y9,) at the scaled proper time 



variables: 



is, to leading order in the perturbation 



^ ^[f{a)]-''\l + 5i5f,59,5ip). 



(3.25) 



We now carry out a gauge transformation which maps 
the shell to a fixed location (note that two different gauge 
transformations are required for the exterior and for the 
interior). For any general gauge transformation = 
— Sf^{x^) we have, to first order in 

where the semicolon represents a covariant derivative 
with respect to the four- metric and gapix'^) is the metric 
tensor in the new coordinate system. We impose that. 



when evaluated at (t, r, 0, (p) = (i, a, 0*, (/3*), the vector 
^" coincides with {St, 6r, 66, 6(p), so that in the new co- 
ordinate system wc will have 



r = a , 



(3.27) 



P* 



where wc are ignoring second-order corrections. The full 
metric in the new coordinate system is 



ga/s 



+ ^KWgaf3 + Sga/3 : 



(3.28) 



where g^^^^ is the static gravastar background metric, 
given by Eq. (|2.ip . The explicit form of and the corre- 
sponding changes in the metric components, Eq. (|3.26p . 
are given in Appendix [XI where the equations of motion, 
as well as the gauge transformation, are presented sys- 
tematically in a multipolc expansion. We then match 
components of gap along the shell, which now simply sits 
at (t, a, 0*, </?*), and apply Israel's junction conditions to 
the extrinsic curvature given by gap- For axial perturba- 
tions these matching conditions read 



^/hhi = 



(3.29) 



For thin shell gravastars h[r) is continuous across the 
shell, implying continuity of the Regge- Wheeler function 
^' and its derivative across the shell [cf. Eq. (|3.9|) ]. In 
more general cases where h[r) may have a discontinuity 
across the shell the axial junction conditions (|3.29p show 
that ^' must also be discontinuous. 

The treatment of polar perturbations is more involved 
and it yields the following relations, determining the 
"jump" of the polar metric functions at the shell: 



m = -87: 



[[Hr]]-2f{a)[[H']] 



l&T:^/J{a){l + 2vl)6Y.. 



(3.30) 
(3.31) 

(3.32) 



Note that /(r) is continuous across the shell [cf. 
Eq. (|2.3[) ]. The parameter Vg in the equations above 
depends on the equation of state on the thin shell, 
e = e(E): 



9S 



(3.33) 



s=o 



and it has the dimensions of a velocity. One might naively 
interpret Vs as the speed of sound on the thin shell and 
require < 1, i.e. that the speed of sound cannot exceed 
the speed of light. Furthermore, for a shell of ordinary, 
stable matter we would have > 0. The standard ar- 
gument used to deduce that > does not necessarily 
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hold when one deals with exotic matter (as in the case of 
gravastars and wormholes), so the specification of upper 
and lower bounds on Vs may require a detailed micro- 
physical model of the exotic matter itself [i^, H^. In 
our discussion of gravastar stability we will consider the 
whole range of Vs, but we will primarily focus on the 
range < < 1. The application of the polar junction 
conditions is more involved than the axial case due to 
their complexity, which arises from the fact that polar 
perturbations couple to oscillations of the shell. 

Here we note that even though we have three quantities 
{K, H, Hi) that satisfy a coupled system of first-order 
ODEs in both the interior and exterior, in each region 
there is an algebraic relation relating the three quanti- 
ties. For this reason we only need to impose two inde- 
pendent junction conditions, and Eqs. (|3.30p - (|3.32p pro- 
vide exactly two independent relations among K, H and 
Hi (after eliminating ST,). Alternatively, the number of 
junction conditions can be obtained considering that all 
metric perturbations can be expressed in terms of \E''" and 
'1'°"*, and that each of these master variables satisfies a 
second-order ODE. More specifically, we use Eqs. (|3.30p - 
(|3.32p to determine two relations among (^''", 9^^'") and 
(^f""', plus the corresponding 5Y,. 



IV. NUMERICAL RESULTS 

Some axial and polar QNM frequencies for a static 
thin-shell gravastar, as computed by the continued frac- 
tion method, are plotted in Figs. [T] and [2l The C-|--f code 
used for the calculations is an update of the Fortran code 
used in Ref. (HO] to verify and extend results on stellar 
oscillations by Kokkotas [51[ and Leins, NoUert and Sof- 
fel [131. For axial modes, our numerical results are in 
excellent agreement with the thin-shell limit of the QNM 
frequencies computed by Chirenti and RezzoUa [36j and 
with Fiziev's calculation of the axial QNMs of ultracom- 
pact objects with a totally reflecting surface (compare 
figures 3 and 4 of Ref. ^). 

To find the QNM frequencies we adopt the following 
numerical procedure. We usually fix = 0.4 and (for 
polar perturbations) we choose a constant value of Ug. 
In the calculations leading to Figs. [T] and [2] we chose, 
somewhat arbitrarily, = 0.1. Later in this section we 
will discuss the dependence of polar modes on Vs- 

As explained in Appendix[Bl within the continued frac- 
tion method the complex QNM frequencies can be deter- 
mined as the roots of any of the n equations fn{^) = 
[cf. Eq. (|B19[) ] , where n is the "inversion index" of the 
continued fraction. To locate QNMs we first fix a value 
of jjL (usually ^ = 0.4). We compute the real and imag- 
inary parts of fn{uj) for a given inversion index n on a 



suitably chosen grid of (w^jjO;/) values, and make con- 
tour plots of the curves along which the two functions 
are zero. The intersections of these curves are used as 
initial guesses for the quasinormal frequencies; more pre- 
cise values are then obtained using MuUer's method [53 |. 
For fixed n (say, n ~ 0) this method singles out some spu- 
rious roots besides the physical QNM frequencies. The 
spurious roots can easily be ruled out, since they are not 
present for different values of n. Looking for roots with 
n = is usually sufficient, but sometimes we get more 
stable numerical solutions for n = 1 and n = 2 when the 
QNMs have large imaginary part (w/ > 1.5 or so). 

The QNM spectrum with [i = 0.4 corresponds to the 
empty (polar) and filled (axial) squares in the left panel 
of Fig. [1] respectively. Starting from [i = 0.4, we follow 
each QNM as ^ and as ^ 1/2 to produce the 
tracks displayed in the figure. For this value of Vs some 
of these tracks end at the origin in the limit /x — > 0, while 
others hit the pure-imaginary axis at some finite limiting 
compactness ^imag (e.g. ^imag — 0.24 for the first po- 
lar mode). The mode frequencies usually move clockwise 
in the complex plane (with the exception of QNMs dis- 
playing "loops" ) as /X is increased. The imaginary part 
of both axial modes (continuous lines) and polar modes 
(dashed lines) becomes very small as ^ — > 1/2, i.e. when 
the gravastar most closely approximates a BH. The be- 
havior is perhaps clearer from Fig. [21 where we separately 
show the real and imaginary parts as functions of /x. 

For both axial and polar spectra the dependence of the 
mode frequencies on the gravastar compactness resembles 
that of "ordinary" ultracompact stars: see e.g. Fig. 3 
of Ref. [5^. Intuitive models that capture most of the 
physics of this problem have been presented in Refs. 
[55[. In their terminology, modes that emerge from the 
origin in Fig. [T] when /x ~ are w-modes or curvature 
modes, roughly corresponding to waves trapped inside 
the star. Modes emerging from the imaginary axis at 
some generally nonvanishing compactness are wn-modcs 
[4?! or interface modes, and they are similar in nature to 
acoustic waves scattered off a hard sphere 1551 The 
only qualitative difference with Fig. 3 of Ref. [53] are the 
"loops" appearing for higher-order w-modes, for which 
we have no analytical understanding. 

The fact that ly-modes are effectively waves "trapped 
inside the star" , while ixin-nrodes are interface modes, 
similar to acoustic waves scattered off a hard sphere, is 
also clear from the behavior of their wavefunctions in 
the stellar interior. The real and imaginary parts of the 
wavefunctions are shown in Fig. [31 (see Ref. [43| for com- 
parison with the wavefunctions of ordinary stars). This 
figure shows eigenfunctions computed at the polar QNM 
frequencies for the first four xw-modes and xfn-modes with 
I ~ 2 and = 0.1. The plot shows that xxi-modes can be 
thought of as standing waves inside the gravastar, and 
that the overtone number corresponds to the number of 
nodes in the real and imaginary parts of the eigcnfunc- 
tion. The situation is different for uiji-modes, where the 
wavefunction has a maximum close to the shell, as ex- 
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FIG. 1: First few axial (continuous lines) and polar (dashed lines) QNMs of a thin-shell gravastar with v1 — 0.1. In the left 
panel we follow modes with i = 2 as the compactness /i varies. In the right panel we do the same for modes with / = 3. Along 
each track we mark by different symbols (as indicated in the legend) the points where /i = 0.1,0.2,0.3,0.4 and 0.49. Our 
numerical method becomes less reliable when 2Mlji is large and when the modes approach the pure-imaginary axis. Numbers 
next to the polar and axial modes refer to the overtone index A'' [N = 1 being the fundamental mode). 




FIG. 2: Real (left) and imaginary (right) part of polar and axial QNMs with / = 2 as functions of /i for v1 = 0.1. Linestyles 
are the same as in Fig. [T] Numbers refer to the overtone index. 



pected for interface modes. 

One of our most important conclusions is that neither 
axial nor polar modes of a gravastar reduce to the QNMs 
of a Schwarzschild BH when fj, — > 1/2. In this limit, 
the real part of most modes is extremely small (much 
smaller than the Schwarzschild result, 2Mu)b. — 0.74734 
for the fundamental mode with Z = 2 [1]). Indeed, the 
QNM spectrum is drastically different from the QNM 
spectrum of a Schwarzschild BH: when yu 1/2 the en- 
tire spectrum seems to collapse towards the origin. This 
is in sharp contrast with the Schwarzschild BH case and, 
as first noted in Ref . [36t , it can be used to discern a very 
compact gravastar from a BH. In Appendix [Cl we prove 
analytically that the QNM frequencies of a gravastar do 
not reduce to those of a Schwarzschild BH as ^ — > 1/2. 



The proof is based on the observation that the Zerilli 
wavefunction for polar modes is continuous in this limit. 

It is clear from the figures that axial and polar modes 
do not have the same spectra for general values of ^. 
However, figures [1] [2] and [4] provide evidence that axial 
and polar QNMs do become isospectral when the gravas- 
tar compactness approaches that of a Schwarzschild BH 
(/i 0.5). An analytic proof of isospectrality in the 
high-compactness limit is given in Appendix [Cl where 
we show that in this limit both the Zerilli and Regge- 
Whceler functions are continuous at the shell. 



From the matching conditions ((X^ . (P^D)) . ([QT]) 
and (|3.32p it is quite clear that polar QNMs (unlike ax- 
ial QNMs) should depend on u^, i.e. [by Eq. (|3.33p ] on 
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FIG. 3: Top row: real and imaginary part of the wavefunction in the interior for the first four w-modes. Bottom row: real and 
imaginary part of the wavefunction in the interior for the first four wn-modes. In both cases we consider polar QNMs with 
; = 2 and = 0.1. 



the equation of state on the shell. This is a new fea- 
ture that does not arise in the case of axial perturbations 
[36l |. The situation closely parallels the ordinary stellar 
perturbation problem [i^, [53|. The role played by the 
equation of state in the dynamical stability of gravas- 
tars against spherically symmetric perturbations was dis- 
cussed in Rcf. [3^. Our calculations extend the consid- 
erations of that paper to nonradial oscillations. 

The Ws-dependence of the modes is studied in Figs.|ll[5] 
andO In Fig. 2] we show the tracks described in the com- 
plex plane by the fundamental polar and axial w-mode 
as we vary the compactness parameter /z. The funda- 
mental axial mode docs not depend on the equation of 
state parameter, as expected, but the polar modes do 
change as a function of Vg. The standard argument used 
to deduce that the speed of sound > does not nec- 
essarily hold when one deals with exotic matter (as in 
the case of gravastars and wormholes). Therefore, for 
completeness, in the left panel of Fig. [4] we compute po- 
lar QNMs when v'^ < 0. Different linestyles correspond 
to different values of u^, as indicated in the legend. The 
solid black line reproduces the fundamental axial lu-modc 



of Figs. [T] and [H The dashed black line corresponds 
to the fundamental polar w-mode for a shell with low 
sound speed (w^ = 0.1), corresponding to the fundamen- 
tal polar w-mode of Figs. [1] and [51 The dash-dash-dotted 
(red), dash-dotted (blue) and dotted (green) lines rep- 
resent a marginally subluminal, imaginary sound speed 
(v^ = —0.99) and superluminal sound speeds (w^ = —2 
and — —4, respectively). Nothing particularly strik- 
ing happens in this regime: QNM frequencies for polar 
and axial perturbations are different in all cases, but for 
large compactness the results become Ws-independent and 
modes of different parity become approximately isospec- 
tral, in agreement with the analytical results of Appendix 
O Furthermore, as \vs\ — > 00 the polar modes seem to 
approach the axial mode. We can perhaps understand 
this behavior if wc think that the shell is effectively be- 
coming so stiff that matter decouples from the space-time 
dynamics, and only the "space-time" character of the os- 
cillations survives. 

The situation is more interesting in the case > 0, 
displayed in the right panel of Fig. U) At first (when 
Vg < 0.5 or so) the modes show a behavior similar to 



10 




FIG. 4: Tracks of the fundamental polar and axial w-modes for different values of the "sound speed" parameter Vs when Vg<0 
(left) and when > (right). Different linestyles correspond to different values of v^, as indicated in the legend. 



that seen for < 0, albeit in the opposite direction (i.e. 
the real and imaginary parts of the QNM frequencies in- 
crease rather than decreasing when \vs \ increases). When 
= 0.7 a cusp develops, and as the speed of sound ap- 
proaches the speed of light (for Vg > 0.9 in the figure) 
the modes "turn around" describing a loop in the com- 
plex plane. The area of this loop in the complex plane 
increases until the sound speed reaches a critical value 
1.013 < w^^.jt < 1.014 (corresponding to Vc-it ^ 1.007). 
For Vs > Wcrit the QNM behavior changes quite drasti- 
cally. The complex mode frequencies still approach the 
LUi = axis clockwise as — > 0.5. However, as de- 
creases the modes approach the axis luj = very rapidly 
along tracks which are now tangent to the lower branch 
of the fundamental axial mode. 



Even more interestingly, when Vg>0 there is also a 
second family of QNMs with very small imaginary part. 
A plot of the eigenfunctions shows that these modes are 
similar in nature to the wn-modes. The trajectories de- 
scribed in the complex plane by some of these weakly 
damped QNMs are shown in the left panel of Fig. [5l 
For comparison, the right panel of Fig. [5] shows some of 
the "ordinary" modes. These ordinary modes are the 
same as in the right panel of Fig. [4l except that this 
time we use a logarithmic scale for the imaginary part. 
The second family of QNMs plotted in the left panel 
has very long damping, and in this sense it is similar to 
the s-modes of "ordinary" ultra-compact stars discussed 
by Chandrasekhar and Ferrari [s^. There is, however, 
an important difference: unlike the s-modes, which ex- 
ist only when a star is extremely compact, the weakly 
damped modes of a thin-shell gravastar only exist for 
small compactness < /icrit- 



In the left (right) panel of Fig. [S] we show the real 
(imaginary) parts of both families of QNMs as functions 
of jj, for selected values of Vg. Both the real and imaginary 
part of the weakly damped modes tend to zero at some 
finite, Us-dependent compactness fi. The range where 
weakly damped modes exist increases with w^. The fact 
that both the mode frequency and its damping tend to 
zero at the critical compactness /icrit suggests that the 
mode somehow "disappears" there, rather than undergo- 
ing a nonradial instability, but this conjecture deserves 
a more careful analytical study. Plots similar to those 
shown in Fig. [S] show that the imaginary part of "or- 
dinary" modes with > 0.84 (vg > 0.92) rapidly ap- 
proaches zero at some finite compactness fi while the real 
part of the modes stays finite. We are unable to track 
QNMs numerically when 2Mu>i < 10~^ using the con- 
tinued fraction method, and in any case we cannot really 
prove by numerical methods that ^ at some finite 
compactness /x > 0. When v1 = 1.0134 the two family of 
modes approach each other along a cusp, but their fre- 
quencies and damping times never cross. The "ordinary" 
family of modes exists all the way up to /x = 0.5 but it 
becomes unstable (in the sense that the imaginary part 
of the mode crosses the real axis with the real part re- 
maining finite) at some finite, Ws-dependent value of the 
compactness /i. 

Summarizing, our numerical results suggest that (1) 
weakly damped QNMs only exist when their compact- 
ness is smaller than some iis-dependent critical thresh- 
old and (2) when > and v1 > 0.84, the imaginary 
part of a QNM crosses the real axis at another critical 
threshold whose value can be estimated by extrapola- 
tion. Nonrotating thin-shell gravastars should be unsta- 
ble against nonradial perturbations when their compact- 
ness is smaller than this critical value. 
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FIG. 5: Left: spectrum of the weakly damped family of QNMs. The vertical line corresponds to twice the orbital frequency 
of a particle in circular orbit at the ISCO: as we will discuss in a follow-up paper, only QNMs to the left of the line can be 
excited by a compact object inspiralling into the gravastar along quasi-circular orbits. In the case Vs — 0.8 the mode "turns 
around" describing a loop in the complex plane. For Vs < 0.8 the modes move clockwise in the complex plane as fi increases. 
For Vs > 0.8 they move counterclockwise and they cross the real axis at finite compactness. To facilitate comparison, in the 
right panel we show again the right panel of Fig. fusing a logarithmic scale for the imaginary part. 




FIG. 6; Real (left) and imaginary parts (right) of the fundamental polar w-mode for different values of the equation of state 
parameter Vs- Diff'erent linestyles correspond to diff'erent values of v^, as indicated in the legend. The horizontal line in the left 
panel corresponds to twice the orbital frequency of a particle in circular orbit at the ISCO: only modes below the line can be 
excited during a quasicircular inspiral. 



The two thresholds are plotted in Fig. [T] Weakly 
damped QNMs with I = 2 exist only below the dashed 
line and, according to our extrapolation of numerical re- 
sults, thin-shell gravastars should be unstable to nonra- 
dial perturbations with I = 2 below the solid line (we 
verified that the instability condition for Z = 3 is less 
stringent than for 1 = 2). The dashed line extends up 
to v'^ < 1.0134, where our numerical search for weakly 
damped QNMs becomes impractical (the modes trace 
smaller and smaller loops in the complex plane and they 
seem to disappear when the compactness is still smaller 
than n = 0.5). 



To validate results from the continued fraction method 
we used an independent numerical approach: the reso- 
nance method [13, [H, , which is applicable to QNMs 
with uji <IC LOfi. The resonance method was first used by 
Chandrasekhar and Ferrari in their analysis of gravita- 
tional wave scattering by constant-density, ultra-compact 
stars [13, HI]. Chandrasekhar and Ferrari showed that 
the radial potential describing odd-parity perturbations 
of these stars displays a local minimum as well as a max- 
imum when the stellar compactness ^ > 0.39. If this 
minimum is sufficiently deep, quasi-stationary, "trapped" 
states can exist: gravitational waves can only leak out 
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FIG. 7: Significant thresholds in the (/^, u^) plane. Thin-shell 
gravastars should be unstable with respect to nonradial per- 
turbations below the solid line, corresponding to "ordinary" 
QNM frequencies whose imaginary part crosses the zero axis 
while their real part stays finite. The dashed line corresponds 
to the "vanishing point" of weakly damped QNM frequencies, 
i.e. to the point where both their real and imaginary parts 
have a zero crossing. We could not find weakly damped modes 
in the region above the dashed line. 
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FIG. 8: The potential governing axial and polar perturbations 
for different values of the gravastar compactness \i = M/a, 
where a is the location of the shell (see [531 1 showing a similar 
plot for axial perturbations of constant-density stars). The 
potential develops a minimum when /i > 0.43. Note that the 
polar and axial perturbations in the interior are both governed 
by the same potential, given in Eg. l3.5l below. 



V. CONCLUSIONS AND OUTLOOK 



to infinity by "tunneling" through the potential barrier. 
The damping time of these modes is very long, so they 
were dubbed "slowly damped" modes (or s-modes) [STj . 

A straightforward analysis of Eq. (|3.5p and an inspec- 
tion of Fig. [5] show that the axial potential for a gravas- 
tar develops a minimum when fj, > 0.43. The compact- 
ness of ordinary stars is limited by the Buchdal limit 
(/i < 4/9 ~ 0.4444), but since gravastars can be consid- 
erably more compact than this limit, s-modes can exist 
all the way down to the "Schwarzschild limit" (/i 0.5). 
These modes can be computed via the continued frac- 
tion method, but since they are long-lived the resonance 
method, which is computationally very simple, provides 
very accurate estimates of their frequencies and damp- 
ing times. We find that the resonance method and the 
continued fraction method are in very good agreement 
whenever the resonance method is applicable. In the 
limit /i ^ 1/2 we have Im(a;) <C Re(a;), and all QNM 
frequencies can be interpreted as "trapped" states. 

The resonance method essentially confirms our contin- 
ued fraction results for modes with ujj <^ oju. In par- 
ticular, it provides additional numerical evidence for the 
conjectured nonradial instability of thin-shell gravastars 
with low compactness. Despite the numerical evidence, 
an analytic confirmation of our estimates of the instabil- 
ity threshold would be highly desirable. 



In this work we have studied the nonradial perturba- 
tions on nonrotating, thin-shell gravastars. It should not 
be too hard to extend our formalism to the more com- 
plex case of five-layer gravastars of the type originally 
proposed by Mazur and Mottola (see [s^ for a treatment 
of the axial case). 

A presumably less trivial extension concerns rotating 
gravastars. Slowly rotating gravastars may be unstable 
against scalar perturbations because of an exponential 
growth of the perturbations due to superradiance, the 
so-called "ergoregion instability" [gO, [6l| . An extension 
of the present formalism can be used to study nonradial 
gravitational perturbations of slowly rotating gravastars 
and to discuss their ergoregion instability, which is be- 
lieved to be stronger for gravitational perturbations [iOl . 
For gravastars this instability is due to superradiant grav- 
itational wave scattering in the ergoregion, so it is essen- 
tially the same as the "w-mode instability" discussed by 
Kokkotas et al. for ultracompact stars [62| . The main 
difference is that gravastars can be more compact than 
constant-density stars, so we may expect the instability 
to be stronger. 

Finally, our formalism can be applied to explore non- 
radial gravitational perturbations of nonrotating worm- 
holes, where the position of the throat plays a role similar 
to the thin shell of a gravastar. Such an analysis could 
confirm or disprove some conjectures on the similarity of 
the QNM spectra of wormholes and BHs [6^ . 

In follow-up work we will extend our study to ex- 
plore QNM excitation by compact objects in closed orbits 
around thin-shell gravastars (see [63.l65j for similar stud- 
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ies considering the excitation of axial modes by scattering 
orbits). According to preliminary estimates by Norte, fol- 
lowing the weak-gravity expansion of Poisson and Sasaki 
[66t (extended by Li and Lovelace 17 1 to general bound- 
ary conditions), the gravitational wave luminosity can 
change dramatically at the resonances, while being very 
close to the BH value in nonresonant conditions. A weak- 
field study should provide a reasonably accurate estimate 
of the orbital frequency at which the resonance occurs, 
but it can only provide a leading-order estimate of the ra- 
diated energy. A more precise characterization requires 
the numerical integration of the perturbation equations 
[67l . m, [1^, [t^ . In a follow-up paper we will compare 
gravitational radiation from EMRIs in thin-shell gravas- 
tar space-times to the BH results of Refs. [7l[ll0,[ll- 
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APPENDIX A: PERTURBATION EQUATIONS 
AND MATCHING CONDITIONS 



In this appendix we develop the formalism to study 
polar and axial nonradial (linear) oscillations of an ob- 
ject consisting of a thin spherical shell separating two 
spherically symmetric regions. Though we are mainly 
interested in thin-shell gravastars, we shall keep the dis- 
cussion as general as possible. We focus on a background 
metric of the form (|2.ip . keeping /(r) and h{r) generic. 

The Regge- Wheeler gauge is incompatible with the re- 
quirement that the shell's world tube sits at a fixed loca- 
tion. We deal with this issue in two steps: (i) we choose a 
coordinate system (system A) such that perturbations in 
both the interior and exterior are in the Regge- Wheeler 
gauge; (ii) we write down the equations of motion of mass 
elements on the shell in this coordinate system (sepa- 
rately for the interior and the exterior), and try to ap- 
ply matching conditions on the boundary of this moving 
shell. We then transform to a new coordinate system 
(system B) in which the shell is fixed, simplifying the 
matching procedure. System B is an auxiliary tool for 
matching: when we consider perturbations of the gravas- 
tar induced (say) by orbiting particles we will mostly use 
system A, which is the usual Regge- Wheeler gauge. 

Step (i) is a straightforward adaptation of formulas in 
Section|TTTl In the Regge- Wheeler gauge, polar perturba- 
tions are defined by three functions K, H and Hi for the 
interior and exterior, respectively. These functions sat- 
isfy an algebraic relation and two coupled ODEs, which 



can be used to reduce the problem to a wave equation. 
This means that we only need two conditions relating 
these quantities from the two sides. Similarly, for axial 
perturbations we have Hq and hi , which satisfy a coupled 
wave equation, and again we need only two conditions 
connecting the interior and exterior. 

Now we go directly to step (ii), and parametrize the 
shell position as in Eqs. p.24p . Let us define 



F = (dF/di) = ^J{a){dF/dT) . 



(Al) 
(A2) 



Note that i is simply a re-scaling of proper time of the 
mass element, and that /(a) is common to the interior 
and exterior, due to the requirement that the induced 
metric is continuous. In the absence of perturbations, t 
coincides with t. Henceforth, we will use (t, 0*,(/J*), to 
parametrize the mass element. As a consequence, the 
four- velocity of the mass element (0* , (/9* ) would be as 



in Eq. (|3.25p . Imposing that QapW 
have to require that 



-1, we actually 



(1 + 25i)gtt{i +St,a + 6r, 6^ + 60, ip^ + d(p) 
W) 

which, to leading order, is 



(A3) 



- 6^wgtt{i,0*,V*) =0. (A4) 

Here S^mgu is the tt-component of the metric perturba- 
tion in the Regge- Wheeler gauge. 

We will now carry out a gauge transformation, both 
in the exterior and in the interior, which maps the shell 
to a fixed sphere. As explained in the main text, for any 
general gauge transformation = a;" — ^"(x'^) we have 

^ {^) (15) 

where we have ignored terms of second order in Not- 
ing that 



(A6) 



we obtain Eq. ()3.26p . This is the desired form, because 
we want to use to transform to a coordinate system 
i'^ where mass elements on the shell arc fixed in spatial 
location and move uniformly in the time direction: i.e., 
for mass elements on the shell, x^^ = (i, a, 0*, (^*). Then 
Israel's junction conditions will be apphed to the Regge- 
Wheeler metric evaluated at a fixed coordinate location 
(i, a, 0*, (y9*), and in terms of the transformation genera- 
tors (which are related to the motion of the shell in the 
Regge- Wheeler gauge). 
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Let us make four consecutive transformations 



e^i) = [0,hir)z{t)Yi„,,0,0], 

?(2) 



sin^ ( 



>(3) 







or, by lowering indices, 

C(0)a - 0,0,0], 

- [0,z(t)ll,„,0,0], 
C(2)q = [0,0, u'(t)y;™|e,u;(t)y;,„|^] , 
C(3)q = 0,0, a;(t)ee^Fi,„l'^,2:(t)e^eyi„J'' 



where the covariant derivative | is defined with respect 
to the 2-dimensional metric 



Gab = de^ + sin^ edip^ 



(A7) 



We impose that, when evaluated at {t,r,6,ip) = 
{t,a,9^,,(p^,), the vector ^" coincides with {St, 6r, 66, Sip), 
so that in the new coordinate system Eqs. (|3.27p will 
be valid. Such a transformation leads to the following 
changes in the metric components: 



2y -jy 



yde yd^ 



Y, 



Im 5 



(A8) 



-fhz 


z 






h' 
h ^ 


zde zd^ 




* 


2rhz 






2rhz sin^ 6 



which are purely polar perturbations; 



Yi^, (A9) 







wYim.e 


wYiyn,^ 


^(2)5/^1/ — 




2 w -y 


— 22!LY, 


* * 















(AlO) 



which are also polar perturbations, and finally 



S (3) 9 til' 





XCQ^Yl^n 
22; Y W 


XC^gYi^j^ ^ 


* * 


X(j)gg 




* * 







(All) 



which are axial perturbations. Here we have defined 

4^mn — -^mlna -^mlma ? (-^-^2) 



where m, n and a run through 9 and (ys. Here eaf, is again 
defined with respect to the metric (|A7p . so that 



(A13) 



In terms of y and the normalization of the 4- velocity 
would be written as 

= -f'{a)h{a)z(t) + f{a)H{i, a) - 2f{a)[~f-\a)y{i)] 
= -f'(a)h{a)z{t) + /(a)F(t, a) + 2y(i) . 

Here we note again that /_)_ = /_ . In this new coordinate 
system the metric is given by Eq. (|3.28p . where 

We are now in a position to match metric components 
along the shell, which simply sits at (i, a, </?*). Of 
course, all of the matching conditions will be expressed 
in terms of the Regge- Wheeler metric perturbations, and 
the motion of the shell in the (internal and external) 
Regge- Wheeler gauges. In the 6 and directions, we 
have 

x+{i) = x-{i), (A14) 
w+{t) = w-{i), (A15) 
2h+ {a)z+ {t) + aK+ (i, a) = 2/i_ (t)z_ (t) + aK^ [i, a) . 

(A16) 

In the to and tip directions we have, in addition, 

V+ (i) = y- {i) , ho+ (i, a) = ho- {i, a) . (A17) 

while in tt direction the matching condition is automat- 
ically satisfied, with 



(A18) 



accurate up to first order in the perturbations. This is a 
consistency check, since we have imposed that the four- 
velocity of mass elements on the shell is 



m" = (i/V7R, 0,0,0) 



(A19) 



which should have a norm of —1. In simplified form, we 
have 

[[x]] = [[w]] = M] = [[ho]] = [[2hz + aK]] 
= [[2ij + fH - ,fhz]] = 0. (A20) 



The symbol "[[ ...]]", as defined by Eq. ([131), gives the 
jump of any given quantity across the shell. 

The four- velocity of mass elements on the shell is given 
by Eq. (jA19p . The surface stress-energy tensor of the 
shell is 

= [j: - e + {Sj: - se)Yi^]ujUk 

- [e + seYi^,]-fjk , (A21) 
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where j and k go through 6 and (p. Now let us try 
to evahiate the extrinsic curvature of the shell at the 
location (i, a,6^,(p^). First of all, we note that 

= (0,1,0, 0)/VF^, 

and the extrinsic curvature is 

where i,j run through t,9,ip. Since Ua only has a 
nonzero r-component the first term vanishes, and we ob- 
tain 



We have 



1 



=r ■ 



From the static configuration of the star, it is easy to 
obtain that 



\K. 



jk\ 



_2/ 



-a^Osin^l 



1 + ^ 



1 + £/l 

^ ^ 2/ 



sin^ Q 



where j , fc = t ^6 ^Lp. From this we have relations (|2.9p 
for static gravastars. 

Now let us focus on first-order quantities. The tensors 
Sjk and Kjk can each be decomposed into six terms: for 
example, 



Sjk 



5(1) 



5(3) 



5(5) 



5(2) sin^ e 



Ylr 



CSC 9dip — sin 



CSC 9d^ 
— sin 



Y\ee Y\g^ 
Y\ip0 Y^^^ 



(A22) 



5(4 
5(5 
5(6 



and 

Ki2) 



fSH + [(9 - 2S)(/F - /iz/' + 2y)] , (A23) 
-a(ais:e + 2hzQ + 0(56) , (A24) 
~^{y + w), (A25) 
I](/io - i) , (A26) 
-29^, (A27) 
-Qx , (A28) 



{H - aK') + fVhz— 2 , (A29) 



a 

a^Vh 



2{Hi + 
/ 

af \ 2K - H 



K' - H' 



^ ■ 2/ y a 
2{h-\L) 2fh + fh' 



a2 af 
ffh' -2{fYh + 2hfr 



2/2 



7^(3) 
^(4) 

^(5) 
^(6) 



a 
Vh 

2Vh 

a 

Vh 



aHi 
— + 

X — 2hQ 
a 

az 



2y + az 
1 



2/ 

+ w 1 + — 
2 V 2/ 



of 
^ 2/ 

/if) - hi 



2 



af 
2/ 



, (A30) 
(A31) 
(A32) 

(A33) 
(A34) 



where Xl = -1(1 + 1). Here A'^'^) an_d K^^^i are axial 
quantities. The junction condition on K^^'^ yield 



\/hhi 



0, 



(A35) 



which together with [[/iq]] = completes the junction 
conditions for axial perturbations (the junction condition 
on (^■' yields an equation of motion for the variable a;) . 
For polar quantities, from [[A'(^^]] = SttS'^^-' we have 



[Whz]] = . 



(A36) 



Then matching [[iir^^']] = STr^^^-* gives us an equation of 
motion for w, while matching [[^"(i'^^]] = SttS^^'"^^ yields 



v-h(^ 



\ a 



2h h' 
a 



Vhz 



(A37) 



Vh (^K' - H' 
Vhd+^M^ 



2Hi 



2f J a 
2h f'h fh' 



Vhz = -167r(59.(A38) 
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The remaining equations are 

[[A']] = -2[[Vh]]Vhz/a:=:8T:T.Vhz, 
fVh 



[[H]] 
SQ 



f 



Vhz 



(A39) 

87r(E - 2e)Vhz , (A40) 
(A41) 



where v is defined as in Eq. (|3.33[) . 

The formalism described above is more general than 
we need for a static thin-shell gravastar. It can easily be 
adapted to more general horizonless space-times and to 
static wormholes. For the Mazur-MottoUa gravastar, we 
have: 



E = 0, e 



[[/'] 



6M 



[[/]] = o = [[/"]], [[/']]- — , [[r 



12M2 



Using the equations of Sec. IIII Al above together with 
these junction conditions we obtain continuity conditions 
for the shell position, [[x]] = [[y]] = [[w]] = [[z]] = 0, and 
the matching conditions for the axial and polar pertur- 
bations functions presented in the main text [Eqs. (j3.29p . 

n^m . <um and 



APPENDIX B: THE CONTINUED FRACTION 
METHOD 

Our numerical search for the QNMs of gravastars is 
based on the continued fraction method, as modified in 
[13, [131 • The QNMs of an oscillating gravastar are so- 
lutions of Eqs. (|3.7|) and (|3.13p satisfying the boundary 
conditions imposed by physical requirements: should 
be regular at the origin, have the behavior of a purely 
outgoing wave at infinity and satisfy the junction condi- 
tions discussed in Section IIII CI The QNM frequencies 
are the (complex) frequencies cu ~ luh + Ilor for which 
these requirements are satisfied. 

The numerical determination of the QNM frequencies 
is nontrivial, especially for modes with large imaginary 
parts (strongly damped modes). The reason is simple 
to understand. Solutions of Eq. (|3.13p representing out- 
going and ingoing waves at infinity have the asymptotic 
behavior - e''*/^ and - e"*^*/^ as oo, 

where r = l/w/ is the damping time. Therefore, identi- 
fying by numerical integration the purely outgoing solu- 
tions (that is, those solutions for which ^f'" is zero) be- 
comes increasingly difficult as the damping of the mode 
increases. The same problem occurs also in the case of 
QNMs of BHs, and was solved by Leaver [t^I. Leaver 
found a continued fraction relation that can be regarded 
as an implicit equation which identifies the quasinormal 
frequencies, thus circumventing the need to perform an 
integration out to large values of r*. This method was 
subsequently adapted to the polar and axial oscillations 
of a star [43, H^. The Regge- Wheeler equation, which 



describes the perturbed space-time outside the gravas- 
tar, becomes 



(Bl) 



with 



and the tortoise coordinate = r -I- 2Mln(r/2M — 1). 
We shall now write the solution of the Regge- Wheeler 
equation in a power-series form as follows. Defining z = 
1 — R2/1', where r = i?2 is some point outside the shell 
of the gravastar, and introducing a function 0(2:), related 
to *(r) by: 



*(r) = (r - 2M)' 



-i2Muj ^ — iLjr , 



(z) = x(r)0(z) , (B3) 



one finds that satisfies the differential equation: 

(co -I- ciz + Ciz^ + csz^)— ^ + (do + d\z + diz^)'^ 



dz^ 



+ (eo + eiz)(j) = . 



dz 
(B4) 



The constants depend only on cu, I and R2 through the 
relations: 

2M 6M 6M 2M 

Co = l-—,Ci = —-2,C2 = l-—,C3 = —, 

U2 -ti2 na -K2 

6A/ „ , „ / 6M\ , 6M 

do = -2iuR2 + — - 2,di = 21-— ,d2 = — , 
-K2 V -K2 / U2 

6M 6M 

eo ^ 1(1 + 1), ei = — — . 

U2 112 

Let us now perform a power-series expansion of (t>{z) : 

00 

0(z) ^ • (B^) 

n=0 

By substituting this expression in Eq. (jB4p . the expan- 
sion coefficients a„ are found to satisfy a four-term re- 
currence relation of the form: 

"102 -f /^ifli + 7iao = , 71 = 1, (B6) 
a„a„+i -I- /3„a„ + 7„a„_i -I- (5„a„_2 = , n > 2 , 

where: 

an = n{n + l)co , n > 1 , (B7) 
f3n = (n — l)nci -I- ndo , n > 1 , 
7n = (n- 2) (71 - l)c2 + (71 - l)di + eo , n> 1 , 
(5„ = (n - 3)(77, - 2)c3 + (77, - 2)^2 + ei , ti > 2 . 

The coefficient ag is a normalization constant, and it is 
irrelevant from the point of view of imposing outgoing- 
wave boundary conditions. The ratio ai/ag can simply 
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be determined by imposing the continuity of ^ and ^' 
at r = i?2, since from Eq. (|B3[) it follows that: 



flO = 


^\z=0 = 




R2 


flO 


*(i?2) 



x{R2) 

*'(i?2) 



R2 - 2M 



*(i?2 



(B8) 
(B9) 



In the axial case, the values of '^{R^) and 5''(i?2) can 
be obtained by the taking the interior solution (|3.8p at 
r = a_ and applying the junction conditions (|3.29p to de- 
termine the wavefunction in the exterior, i.e., at r = a+. 
From then onwards, we can numerically integrate the 
Regge- Wheeler equation (|3.13p up to r = R2. The re- 
maining coefficients can then be determined by recur- 
sion from Eq. (jB6p . In the polar case we proceed in a 
similar way: we obtain the Zerilli function and its 
derivative at r = a_|- by imposing the matching condi- 
tions ((3?30l) - ([3?32|) . Then we use Eq. (|3?24l) to obtain the 
corresponding Regge- Wheeler function at r = a_|_, inte- 
grate forwards to find ao and ai, and finally obtain the 
remaining coefficients by recursion. 

To apply the continued fraction technique, it is easier 
to consider three-term recurrence relations. Leaver has 
shown that the four-term recurrence relation (|B6P can be 
reduced to a three-term recurrence relation by a gaussian 
elimination step [t^. Define: 



f3o 



ao 



do = -1 



(BIO) 



where ai/ao is obtained numerically from Eq. (|B9[) . Now 
set: 

an = an, Pn = (3n , = 0, 1) , 

> = 7n , {n=l), (Bll) 



and for n> 2: 



/3n=/9n 



Oin-lSn 
7n-l 



Pn-lSn c n 

7n = 7« : , On = . 

7ii-l 



(B12) 



By this gaussian elimination, Eq. (jB6|l reduces to: 

+ Pnttn + %an-l =^ . (B13) 

The elimination step is not as trivial as it may seem, 
because in the process one of the three independent so- 
lutions to Eq. (|B6p is lost. It can be shown that this 
solution is not relevant for our purposes [i^ . 

We now turn to investigating the asymptotic behavior 
of the coefficients a„ in the expansion (|B5|) . Let us make 
the ansatz: 



V fln+i h k 

hm = 1 H -J- H — 



(B14) 



Dividing Eq. (|B6[) by n?an, keeping terms up to ^ n~^l'^ 
and equating to zero the various terms in the expansion 
in powers of ?7~^/^ we find the relations: 



(B15) 



Co + Ci + C2 + C3 = , 

2co + ci - C3 = , 

= HujRi , 
3 

A: = - - + 7tj(i?2 + 2M) . 



The first two of these equations are identities. Substitut- 
ing the second pair of equations in Eq. (|B14p we get: 



lim a„ 

n — ^oo 



-3/4+iw(fl2+2J\f)g±2\/2iw_R2ii 



(B16) 



According to a definition given by Gautschi [73| , the so- 
lution of Eq. (jB14p corresponding to the plus sign in 
Eq. (|B16P is said to be dominant, whereas that corre- 
sponding to the minus sign is said to be minimal [tJ. If 
we select the minimal solution the expansion (|B5p is ab- 
solutely and uniformly convergent outside the star, pro- 
vided that we choose R2 such that R2/2 < a < R2 and 
i?2 > 2. Furthermore, according to Eq. (|B3p . the solution 
to Eq. (|Bip behaves as a pure outgoing wave at infinity, 
i.e. it is the a QNM wavefunction. Thus, the key point is 
to identify the minimal solutions of Eq. (|B14p . Accord- 
ing to a theorem due to Pincherle [tJ , if Eq. (jB14p has 
a minimal solution then the following continued fraction 
relation holds: 



«]_ _ -71 ai72 '3273 

«o Pi- $2" h- 



(B17) 



where the continued fraction on the RHS is convergent 
and completely determined since the coefficients d„, /3„ 
and 7„, defined in eqs. (|Blip . (|B12p are known functions 
of u). Moreover, from Eqs. (|B9p and (jBlOp it is apparent 
that the dependence on the stellar model is all contained 
in the ratio ai/ao. Keeping in mind the definitions (jBlOp . 
Eq. (jBlTp can be recast in the form: 



= fo{uj) = Po 



ao7i q;i72 ^273 
Pi- P2- Ps- 



(B18) 



Using the inversion properties of continued fractions [78| , 
the latter equation can be inverted n times to yield: 



o = /„H = p 



ctn-iln a„-27K-i q;o7i 

Pti-l— Pn~2— Po 
Q!„+i7„+2 d„+27„+3 



•n+2- 



P 



(B19) 



'n+3" 



for n = 1, 2, .... These 77 conditions arc analytically equiv- 
alent to Eq. (jBlSp . However, since the functions /n(w) 
have different convergence properties, each of them is 
best suited to find the quasinormal frequencies in a given 
region of the complex uj plane. This is the main reason 
for the accuracy and flexibility of the continued fraction 
technique. 



18 



APPENDIX C: HIGH-COMPACTNESS LIMIT 



terms in a — 2M, 



To investigate the behavior at the surface of the gravas- 
tar in the high-compactness Hmit, we use the z —>■ 1 — z 
transformation law for the hypergeometric function [i^ . 

F(a, b, c, z) 

^ ^ r(a)r(5) 

xF{c~a, c— &, c— a — 6+1,1 — z) 
r{c)r{c - a - b) 



T{c- a)T{c-b) 



F{a,b,-c+a+b+l,l-z). (CI) 



Using Eq. p.Sp in the Hmit when C ^ 1 and r — > a we 
get 



2(a- r) 



T{1 + |)r(-i2Mw) 



l+^-^2Mt£^ 



2 (a - r) 



V{1 + §)r(i2A/w) 



(C2) 



Within our conventions the first term is in-going, while 
the second term is out-going near the surface. So it is 
clear that in this regime both in- and out-going modes are 
present, and QNMs do not reduce to the Schwarzschild 
QNMs (which require only in-going waves) . Furthermore, 
we can clearly see that in-going and out-going waves al- 
ways have the same magnitude: the gravastar appears 
like a reflecting object, as suggested by [t^. Because 
this reflection happens in a polar coordinate system, it 
can simply be interpreted as due to the fact that waves 
going into a lossless gravastar will re-emerge without loss. 
Nevertheless, such a behavior already supports the con- 
clusions of Ref. [sO, [lOl, which showed that (for scalar 
fields) the ergoregion instability is more effective when 
the surface of the compact objects behaves like a "per- 
fect mirror" in this sense. 

It is easy to show that in the high-compactness limit 



(C3) 



a-2M 

Solving for the metric quantities we find, up to dominant 



if(a_ 



l{l + 1) + 2iauj 



2a 

+ 1) + 2uja 
'"^ 2(a - 2M) 
to a(ui a — i) 
a - 271/ 

a — i){4:M + iuja'^) 
(a - 2M)2 
to a{LU a — i) 

a-2M 
io{ui a - 2i){2M + 11^0^) 



In the exterior we get 



(a - 2M) 



*(a_), (C4) 

*(«-), (C5) 

*(«-), (C6) 

vl/(a„), (C7) 

*(«-), (C8) 
*(a_). (C9) 



K{a+) 
K\a+) 

Hi{a+) 
H[{a+) 
Ho{a+) 



*(a_). 



1(1 + 1) + 2iau 
2a 

+ 1) + 4wa 
2{a - 2M) 
bjM{i + AMoj) 

a - 2M 
Muj{2Muj + - AiMio) 

(a - 2M)2 
Muj{i + AMuj) 
a-2M 
H[{a+). 



(CIO) 
(CU) 

*(«-), (C12) 

^(a-), (C13) 

vE'(a_), (CM) 
(C15) 



Notice that, even though Hi is not continuous at r = a, 
the Zerilli function is. Indeed, we get 

Z°"*(a+) = *(a_), ^{a-)^^{a+). (C16) 

Thus, we conclude that in the high-compactness limit, 
the master wavefunction for polar perturbations is con- 
tinuous across the shell. A trivial extension of the known 
Schwarzschild results then shows that polar and axial 
perturbations are isospectral for large compactness, i.e., 
when a 2M and — > 1/2. 
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